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We extract the next-to-next-to-leading order results for spin-flip generalized polarizabilities (GPs) 
of the nucleon from the spin-dependent amplitudes for virtual Compton scattering (VCS) at Oip^) 
in heavy baryon chiral perturbation theory. At this order, no unknown low energy constants enter 
, the theory, allowing us to make absolute predictions for all spin-flip GPs. Furthermore, by using 

' constraint equations between the GPs due to nucleon crossing combined with charge conjugation 

I symmetry of the VCS amplitudes, we get a next-to-next-to-next-to-leading order prediction for one 

, of the GPs. We provide estimates for forthcoming double polarization experiments which allow to 

. access these spin-flip GPs of the nucleon. 

OX). 
< 



On 
O 



PACS numbers: 13.60.Fz, 12.39.Fe, 13.40.-f,14.20.Dh 



I. INTRODUCTION 

. Over the past decade, the virtual Compton scattering (VCS) process on the nucleon, accessed through the ep epj 
reaction, has become a povirerful new tool to study the internal structure of the nucleon both at low and high energies 
(see 0] for a review). At low virtuality and energy, the outgoing real photon plays the role of an applied quasi- 
■ static electromagnetic field, and the VCS process measures the response of the nucleon to this applied field, which 
can be parametrized in nucleon structure quantities, termed generalized polarizabilities (GPs) 2J. In this case, the 
virtuality of the initial photon can be dialed so as to map out the spatial distribution of the electric polarization, the 
magnetization, or the spin densities of the nucleon Because the chiral dynamics plays a dominant role in this 
regime, at low virtuality. Heavy Baryon Chiral Perturbation Theory (HBChPT) provides a natural way to study the 
O i' GPs of the nucleon. 

' ! The spin-dependent VCS amplitudes in HBChPT at order 0{p^) have been calculated in Refs. Hll- The HBChPT 
calculation for the spin-dependent VCS amplitudes has recently been extended to 0(j5*) in Ref. At this order, no 
unknown low-energy constants enter the calculations so that one can extract the next-to-leading order (NLO) results 
• for the spin-flip GPs of the nucleon, without free parameters. In Ref. i^], it was possible to extract the next order 
. 5^ I corrections to three of the seven spin-flip GPs in HBChPT by identifying the terms linear in the final photon energy 
i in a low-energy expansion (LEX) of the spin-dependent VCS amplitudes. 
^ ' In this paper, we demonstrate that by performing higher order expansions of the VCS amplitude with respect to 
the final photon energy, i.e. by including the terms quadratic in the flnal photon energy, one can furthermore extract 
the NNLO results of three more spin-flip GPs and even make a N^LO prediction for one particular CP. In this way 
we are able to obtain the next order corrections to all seven spin-flip GPs in HBChPT, without free parameters. 

On the experimental side, the measurement of the VCS process became only possible in recent years, with the advent 
of high precision electron accelerator facilities, as these experiments involve precision measurements of small cross 
sections. At low photon energy and virtuality, first unpolarized VCS observables have been measured at the MAMI 
accelerator at a virtuality = 0.33 GcV^, and recently at JLab Q at higher virtualities, 1 < < 2 GeV^, and 
data are under analysis at MIT-Bates All those experiments measure two combinations of GPs. To disentangle 
the four independent spin-fiip GPs, requires to perform in addition double polarization experiments. In this paper, 
we show predictions for such double polarization asymmetries and demonstrate their sensitivity to the spin-flip GPs. 

The outline of this paper is as follows. We start in Section ^ by reviewing how the virtual Compton scattering 
amplitude at low photon energies can be expanded in terms of GPs. We subsequently review the status of the 
HBChPT calculations for the spin-flip GPs. In Section IIIII we outline how we can extract four spin-flip GPs at 
next order in HBChPT by performing a quadratic expansion in the final photon energy. In Section llVI we compare 
our predictions for the four independent spin-flip GPs at next order in HBChPT with the leading order (non-zero) 
expressions. Furthermore, we demonstrate the sensitivity of double polarization asymmetries to the spin-flip GPs. 
Finally we give our conclusions in Sectional Technical details on crossing symmetry relations for the VCS amplitudes 
can be found in an Appendix. 
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II. VIRTUAL COMPTON SCATTERING AND GENERALIZED POLARIZABILITIES 



We begin by specifying our notation for the VCS process: 

7*(£i,g) + 7V(pi) -.7(e2,'z') + ^b2)• 



(l) 



The VCS amphtudes are calculated in the cm. frame and we choose the Coulomb gauge |0j. In the VCS process, the 
initial spacelike photon is characterized by its four momentum q = {uj, q), virtuality = —q^, and polarization vector 
ei = (0,£i). The outgoing real photon has four momentum q' = (w' = l^'l,?') and polarization vector 62 — (0,e2)- 
We define q = \q\, and denote as the scattering angle between virtual and real photons, i.e., cos 6' — q ■ q . The 
polarization vector ei of the virtual photon can be decomposed into a longitudinal component sil = (ei ■ q)q and a 
transverse component Sit- For further use, we introduce the virtual photon energy in the limit lu' ^ : 



u;o=uj{u;' = 0,q) = M^- \jMj^ + q^, 

= -f/(2MN)+0{l/Ml), 

with Mn the nucleon mass, and where the last line in Eq. ((21 indicates the heavy baryon expansion. 
The VCS amplitude Mvcs can be expressed in terms of twelve structure functions as (J| : 



(2) 



M 



VCS 



£3 • £it)Ai + (£2* • q){£lT ■ q')A2 + id ■ (£2* X £it)A3 + 1(7 • (q' X q){el ■ £it)^4 



+ id ■ (£2* X q){eiT ■ + iS ■ (£2* X (7')(eiT • 9') A " ' (^it x g')(e2 ' q)M - io ■ (eit x q){e2 ■ q)A^ 

+ (£iL • q) (£2* • q)A^ + ia ■ {q' x q){e2 ■ q)AiQ + iff ■ {e^ x q)Aii + iff ■ (£2* x q')Ai2 } X» ■ (3) 

To extract the GPs, we first calculate the complete fourth order VCS amplitudes Ai in HBChPT, and subsequently 
separate the amplitudes in a Born part Af °™, and a non-Born part Ai as : 



A{u:\ q, 9) - Af q, 9) + A,{J , g, 9), 



(4) 



In the Born process, the virtual photon is absorbed on a nucleon and the intermediate state remains a nucleon. We 
calculate them, following the definition of Ref. , from direct and crossed Born diagrams with the electromagnetic 
vertex given by : 



qv 



2M 



N 



(5) 



where Fi{F2) are the nucleon Dirac (Pauli) form factors respectively. The non-Born terms contain the information 
on the internal structure of the nucleon. In particular, we are interested in this work in the response of the nucleon 
to an applied quasi-static electromagnetic dipole field. This is accessed by performing a low energy expansion (LEX), 
in the outgoing photon energy, of the non-Born VCS amplitudes, and by selecting the term which is linear in the 
outgoing photon energy, given by : 



9A 



H^\q,0) 



= S,{q)+'P^{q) COS 9. 



(6) 



In Eq. ((BJ, the nucleon structure quantities Si and Vi can be expressed in terms of the GPs of the nucleon, as 
introduced in which are functions of q and which are denoted by P*^'' ^ '^^^^{q). In this notation, p (p') refers to 
the electric (E), magnetic (M) or longitudinal (L) nature of the initial (final) photon, L {L') represents the angular 
momentum of the initial (final) photon, and S differentiates between the spin-flip {S — 1) and non spin-flip (S — 0) 
character of the transition at the nucleon side. Restricting oneself to a dipole transition for the final photon (i.e. L' 
= 1), angular momentum and parity conservation leads to 3 scalar and 7 spin GPs 0. The 7 spin GPs are obtained 
as 



a a 

m 

p{MlX2)l 
p(Ll,M2)l 
p(A/l,2)l 



-2V2 1 

3V3 qi^a V Em 
V2 



zll. Mes p(Mi.Lo)i^^± Me. 

3 ujoy En V En 



Si 



2. ■ 



3g2 



N 



5 p(Ml,Ml)l 



N 



2 



N 



N 



[S7 - Si] 
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Mn 
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S3 + t:Ss — Su 



(7) 
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where En = \/ + 9^ ■, and where the GPs denoted by P correspond with mixed electric and longitudinal multipoles, 
as introduced in It has been shown [1(1( 1 that nucleon crossing combined with charge conjugation symmetry of the 
VCS amplitudes provides 3 constraints among the 7 spin GPs : 

Si = 0, ^3 = —St, Sio — Si2 = —Sii, (8) 

leaving 4 independent spin GPs. By expanding ujq as in ((SJ, it is obvious that the relations JSJ connect quantities of 
different order in the heavy baryon expansion. These constraints have been verified in HBChPT in Ref. [g. 

By calculating the VCS amplitudes in HBChPT at 0{p^^), one extracts from Eq. lO the following expressions for 
the GPs at LO (see Ref. [|) : 

^p(Ml.Ml)l(-)]^0 ^ [p(Ll,Ll)l(-)]^0 ^ 

t^'""""-'l"'-i3v^{'-OT""(f)}- 

where w = q/m-^, with TOtt the pion mass. Furthermore, throughout this paper we use the values : qa — 1.267, 
F.^ = 0.0924 GeV, and ^ 0.14 GeV. 

The VCS amplitudes in HBChPT at 0{p^), yield the GPs at NLO, from Eq. (0 (see Ref. 0) : 

rp,Mi..2)i( ^niv^o ^ g\ \ 4 

^ 12V67r2F2Q2 I vjyflPTl V 2 y 

L 12V37r2F^ \ V 2 y / 

rp(Mi,Mi)i.->,iJViO g\ f +4 . , -1 /«; 

24yT07r2F2Afjvg2 \ w^w^ + 4 

, JVLO 



[p(il,il)l(^-)]"-- ^ 
^p(Ll,A/2)l(-)^ 



(Li,M2)i,-x]^io _ gj </ 



96\/27rF|g2 M]v 

< 1 ^ r + ( T)tan — + ^3 ^( ^)tan — ^ 

\2w «;m2+4 M lu^ ^ 12/ L2w M ^ V2yjj 



fll 2w^+4 ,25 9, -1 /«;\ [3 ,5 3, 
\2^' ^(u,2+4) '(T + I^^^^" l2j+"n2^"^4 + I^)*'^^ 



(10) 



One can get more predictions by use of the crossing relations (jH)), which hold in general in a relativistic quantum 
field theory, and which were calculated in HBChPT to ©(p^) in Ref. 0. By plug ging in the fourth order amplitudes 

in JHJ, we can extract fifth order predictions from ^g^'' — —2Mis[ I qS\^\ and ^^q'' — s'^ = —2Mn /qS[f\ In this way, 
these relations allow us to extract two spin GPs at next-to-next-to-leading order (NNLO) : 

rp(Mi,Min( ^nA'iv^o ^ -g\g f 3 ^ 5 6 x M _ 1 2 , | 



L ^^^J 487T^FiMN \ wV^^i;^T4 \2j j ' 

It is surprising that the third order calculation (5| was able to obtain some of these NLO and one NNLO results! 
Actually, the_NLO results for p(Mi,L2)i ^^^^ p(Mi,Lo)i ^^^^ obtained in Refs. 0,13, by performing the LEX of the 
amplitudes to second order in to', and, by isolating two terms at order w'^ which depend on those spin GPs. 
Furthermore, they used the crossing relations ||SJ) to obtain the NLO results for p(*^i'^^^i)i and p(*^i'2)i^ ^^yi as 
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the NNLO result for p(LhLi)\ 

In the present work, we generahze this method and obtain the NNLO result for p(A/i,i2)i^ p(Afi,LO)i^ o^j^j p(Mi,2)i^ 
and even the N^LO result for The first step is to define the following quantities: 

= 2a,(g) + 2/3,(g)cos6l + 27,(g)cos^6'. (12) 

Our goal is to express a^, (3i and 7^ in terms of GPs. To achieve this goal, one needs to start from the covariant 
virtual Compton scattering tensor. The VCS amplitudes are obtained as the contraction of the VCS tensor M^i^ with 
the polarization vectors of the photons, evaluated between the nucleon spinors in the initial and the final states, 

12 

Mvcs = -ie^uijpi) ^ ei^el^pl"" fi{q^,q ■ q, q ■ P)u{p2), (13) 
1=1 

where P = pi + v^j-fi are the VCS amplitudes, and p^" are the gauge-invariant independent tensor structures for 
VCS, as given by [10|. The key point of our method is to connect the covariant amplitudes fi and the c. m. amplitudes 
Ai defined in Eq. (PJ- Subsequently, we will perform a low energy expansion of the amplitudes Ai and sort out the 
relations between fi and the ai,/?i and 7^ defined through Eq. I|12|) . In addition, the GPs can also be expressed in 
terms of the covariant amplitudes fi. Comparing both results will enable us to express the higher order coefficients of 
the LEX in terms of GPs. The advantage of working with the covariant Compton tensors is that symmetry properties 
due to photon crossing as well as nucleon crossing combined with charge conjugation, as derived in Ref. 10] and 
detailed in the Appendix, are immediately manifest. 



d'Ai{cu\q, 



III. CALCULATION OF GENERALIZED POLARIZABILITIES AT NNLO AND N^LO IN HBCHPT 

To extract the GPs at NNLO, we start by connecting both sets of VCS amplitudes Ai and fi of Eqs. Q and H13|l . 
To this goal, one has to expand £ifj.£2uPi'^ ^1- P3|l in lu' in the cm. frame and also expand the non-Born parts of 
the amplitudes fi in uj': 



f, 



non— norn 



Mq) + ^'- [9: {q) + h^{q) cos e] + o{uj'^). 



(14) 



For the new predictions of three NNLO GPs and one N'^LO CP which we will derive, we only need the amplitudes 
^10,^11 and A12. The expressions for these amplitudes including terms up to order w'^ can be found, after some 
algebra, as [20| : 
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(15) 



5 



+ {—2ljq +2LUoqcos9)hii + ( —2Mm'^o + 



Mnujo 
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(16) 
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(17) 



We can simplify these expressions by using the crossing symmetry relations for the amplitudes fi, as given in the 
Appendix, which lead to /4 = /s = fio = 0. It is now straightforward to read off from Eqs. (|15I - I17|I the coefficients 
Si of Eq. ((HJ, which appear in the LEX of the cm. amplitudes Ai, and express them as combinations of the fi : 



Sio — 
Su = 

Si2 = 



2q 



— -/7 - 2ujoqfii : — /i2 
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(18) 



By applying the defining relations of Eq. |7| for the GPs, one then obtains 
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(19) 



In general, the coefficients of order w'^ are combinations of /, gi and hi. On the other hand, the dipole GPs defined 
through Eq. 10 are all combinations of /. Actually gi and hi are related to higher order polarizabilities {L' > 1), and 
they can be defined in a similar way as in Ref . j3| ■ We leave the study of such higher order GPs to a future work and 
concentrate in this work on new predictions for dipole {L' — 1) GPs. 

To compare with the amplitudes obtained from HBGhPT, it is necessary to make the heavy baryon expansions on 
Eqs. ((T5l - [T7|) . Note that not only luq has to be expanded in q/M^ as in Eq. (0), but also /i, gi and hi. Denoting any 
quantity proportional to (I/Ma?)" by a superscript (n), we can express fi as : 



/. - /f ^ + /f ^ + 0{llMl), i = 1, 5, 6, 7, 9, 11. 



(20) 



The amplitude /12 as well as the combination 2/g + /g start with a contribution at n = 1 otherwise the amplitudes 
As and A12 would start at n = — 1, which is obviously impossible. Also /2 starts at n = 1 (otherwise Ag would also 
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start at n = — 1), so that we can write : 



/, -/f^+0(l/M2,), z^2,12. 



(21) 



The coefBcients in the heavy baryon expansion are furthermore constrained by the crossing symmetry relations for 
the amphtudes fi (for details, see the Appendix). We can then separate the leading and subleading terms in the 
heavy baryon expansion for the amplitudes ^lo, ^ii, and A12 as : 
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From Eq. (|22l) one can now identify the coefficients , (3i and 7^ appearing in the LEX of Eq. p2|) as follows : 
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The GPs of Eq. p9|l can now be expressed in terms of the above LEX coefficients. The power counting of the 
GPs needs some word of explanation though. From HBChPT, the leading order amplitudes Si ^ C(l) and 
S^"^^ ^ 0(1/Afjv). Because loq is a quantity of order 0{\/Mj^), one remarks from Eq. {Tj) that the three GPs 
p(Mi,L2)i^ p(Li,Li)i p(Mi,L0)i should Start at 0{Mn). Actually s'f^ = s'f^ = = 0, and consequently these 
three GPs are non-zero only at NLO : 
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To evaluate the GPs piMi,L2)i ^^^^ p(Af i,LO)i nnLO, requires the knowledge of the VCS amplitudes at O(p^) 

In general, the amplitudes at 0{p^) in HBChPT contain 



in HBChPT through the coefficients a%^^ and afa^^. 

two parts. The first one consists of two-loop diagrams and one- loop diagrams with one vertex from the third-order 

(3) 

Lagrangian >C^^. The second one consists of one-loop diagrams with vertices from the heavy baryon expansion. The 
main difference between both contributions is that the former ones contain no 1/Mn factor whereas the latter ones 
contain a l/Af^ pre-factor. We will present here the results due to the one-loop diagrams with vertices from the 
heavy baryon expansion. Taking into account these contributions will precisely restore the symmetry property due 
to nucleon crossing combined with charge conjugation, which relates the VCS amplitudes at different orders in the 
heavy baryon expansion. The new analytic structures which one can expect to arise at 0{p^) due to two- loop effects 
have then to satisfy the crossing relations separately. 

We find for a^^*^ and a^2^'-' by explicit calculation in HBChPT the following results : 



a 



NLO 



= 0, /3(f = 0, 



10 



NLO 
"12 



9\ 



4 6 + 15-u;2 
3 ' 



2ii;4 



+ 7:[(1 + K,„) - (1 + k,)t3 



MMnt^F^ \ \_w Aw + 3 Aw + 

I ^ , , r 8 w ,„ 16 , ^ fw\ 

- + 2+ — tan-M- 

w 6 V 2 / 

f r-7 2w 14 , -I fw\ 

.)tan-(|) 



(12 



9\ 



I, 



-3 
w 



w , 1 

6+^(2 



-1 w ,5 2 , (w 
--77+ x + — tan M- 
D \S 2 \2 



(25) 



where — 3.70 ( Kg = —0.12 ) are the nucleon isovector (isoscalar) anomalous magnetic moments respectively. 
Combining Eqs. (|24|l and H25|) . allows us to obtain the NNLO results for two GPs : 
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A third NNLO prediction follows by applying the crossing symmetry constraint of Eq. 0, i.e. Sio 
which yields the relation between 3 spin GPs : 
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After the heavy baryon expansion, this becomes : 



6Mn [p(^^^'^^')^]'^^0 ^ ^^p{MlX0)l-^NNLO _|_ ^1 :i-2^p{Ml,L2)l^NNLO 



(29) 



from which we obtain : 

N^LO 



p(Ll.Ll)l(-) 



9a? 



15 4-6^2 6 + 15u;2 + 2u)'' 
2u) + 3 • - 4 

W 4:11] + 4:11) + 



.,30 -5,..- (I) 



+ [l + Kv] 
+ [(1 + Ks)t3 



15 w ,3 15 , 1 / 
2w 4 M 7«2 ^ V 2 



27 w , 27 27, 1 1 

h( ^)tan-M— } 

2w 4^4 w^' ^2 /J J 



(30) 



Similarly, by applying the crossing symmetry constraint of Eq. Q, i.e. 54 — 0, which yields the relation between 3 
spin GPs : 



3^^p(Ml,L2)l(-) ^ p(Ml,Afl)l(-) ^ /5 -2p(A/l,2)l ^ 



(31) 



one can get the NNLO result of p(*^i'2)i 

NNLO 



p(Ml,2)l(-) 



96VlOAf|,7ri^2^ 
4 



-5 4-6w2 4 6 + 15u;2 + 2w'^ ,29 12 



w 4w + 



4w + 



\ tan 



+ [(1 + K^) - (1 + Ks)t3 







(1)] 


+ T3 



11 2 , 

+ + — tan- 

w 2 w' 



.(32) 



At the real photon point ( g = 0), we can express the GP p(*^i^^2)i(^Q-j terms of a sum of two spin polarizabilities 
(72+74 ) , introduced by Ragusa ^3] i as : 

^, + ^, = _«_^p(Ml,L2)l(o)^ 

where = e^/(47r) = 1/137. The result of Eq. (|?7|l then allows to extract the result for 72 + 74 at NLO as : 



(72 + 74) 



NLO 



-3 + 2(1 + k„)-2(1 + k,)t3], 



(33) 



384TTF^m^MN 

which exactly reproduces the NLO calculations for 72 + 74 of Refs. 

From our results, it is also worth noting that the crossing relations of Eq. (|SJ| allow to classify the GPs into two 
groups : 

group 1 = 



^ |p(Afl,L2)l p{Ml,LO)l p{Ml,Ml)l p(Afl,2)l| 
grOUp2 = |p(il,il)l,p(il,M2)l p(Ll,l)l| 



It is interesting to observe that their analytical forms alternate from one order to the next, as noticed already in 
Ref. At LO one has : 



At NLO, one has : 



jroup 1 : — tan 

TT 



jroup 1 : sinh 



2m, 



; group 2 : — sinh ^ 



2m, 



2™, 



^1 -1/9 
jroup 2 : — tan 

TT \ 2m, 



At NNLO one again has the same analytical structure as at LO, etc. One therefore observes that the alternating 
analytical forms teach us how the crossing relations ijSJ), connecting the GPs between two different orders in HBChPT, 
work. The same alternating analytical form were observed when comparing leading and next-to-leading order HBChPT 
results for the forward spin polarizabilities in Ref. |l4j . 
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IV. RESULTS AND DISCUSSION 



The main result of this work are the new NNLO predictions of Eqs. (|26l I27() and (|32|1 for the spin GPs p(Ml,LO)l 
p(Mi,L2)i^ and f(*^i^2)i^ ^1^^ jv^Slq prediction of Eq. (HU for the spin GP pi^hLi}!^ 



p(u,u)i P^"'-"'^' 





0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



(GeV^) 



(GeV^) 



p(U,M2)1 ^ ^ Q-3 ^^4>) p(M1.L2)1 ^ ^ q-3 ^^4>) 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



(GeV^) 



(GeV^) 



FIG. 1: Results for the spin GPs in HBChPT. For the GP p(^i'"2)i^ ^^^^ LO (dotted curve) and NLO (dashed curve) 
resuhs are shown. For the GPs p(Mi,Mi)i j^^^j p(Mi,i2)i^ ^^^^ ^-^q resuhs (dotted curves) and the NNLO resuhs 
(dashed curves) are shown. For the GP p(^i.^i)\ the NNLOJclotted curve) and N^LO (dashed curve) results are shown. For 
comparison, we also show the dispersive evaluation of Refs. [islllfil (solid curves). 



In Fig.^ we compare these results for the GPs p(^^i.^2)i ^^^^ p(Lixi)i ^j^j^ ^j^^^-j. i^^ggt oj-jer non-zero expressions, 
and display as well the next order results for the GPs pl^-i.^fa)! p(Mi,Mi)i ^ which were calculated before in 
Ref. 0. Note that the three remaining spin GPs of Eq. iQ can be expressed in terms of these four independent 
spin GPs by making use of the crossing symmetry relations of Eq. ((SJ. By comparing the lowest order with next 
order results in HBChPT for each of the GPs, one observes larg e corrections at next order. We also compare these 
results in Fig. with the phenomenological estimate of Refs. |l5l lT^ using dispersion relations. One sees that for the 
GPs p(i'i.Af2)i^p(j\/i,Mi)i g^^^ p(Li,Li)i ^.jjg corrections to the leading order bring the HBChPT results towards the 
phenomenological estimates, though overshoots it for For the GP p(A^i.^2)i^ ^y^^ sizeable large correction 

at the real photon point brings the HBChPT results close to the dispersive results but shows a much slower fall-off 
with than the dispersive estimate. The comparison in Fig.ndearly indicates that a satisfactory description of spin 
GPs is still a challenging task. To completely disentangle the four spin GPs experimentally, requires four independent 
observables, which we discuss in the following. 

It was shown in Ref. that GPs can be accessed experimentally by measuring the ep ep7 reaction and performing 
a low energy expansion in the outgoing photon energy w'. In particular, the VCS unpolarized squared amplitude J\A°^^ 
takes on the form Q : 



+ M°"'' + 0(c^')- 



(34) 



Due to the low energy theorem (LET) 
determined from the Bethe-Heitler -I- Born (BH 



the threshold coefficients A^f'J^ and M^^f are known and are fully 



B) amplitudes. The information on the GPs is contained in A^q'^ 
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which contains a part originating from the BH+B amphtudes and another one which is a Hnear combination of the 
GPs, with coefficients determined by the kinematics. The unpolarized observable Ai^^^ can be expressed in terms of 
three structure functions Pll{<i), PTxiq), and PLriq) as : 

- -^r^^ = 2i^2|«i [ePLLiq) - Prrm + («2 - j^a) V^^{l + e}PLT{q)^ , (35) 

where K2 is a kinematical factor, e is the virtual photon polarization (in the standard notation used in electron 
scattering), and ui, W2, are kinematical quantities depending on e and q as well as on the cm. polar and azimuthal 
angles {6^^ and (jj, respectively) of the produced real photon (for details see Ref. j3|)- The three unpolarized observables 
of Eq. (|35|l can be expressed in terms of the independent GPs as ;2J • 



Pll = -2V6MjvG£;(Q2)p(ii.ii)o , (36) 



f 

UJQ 

Plt = xll^GEmP^'''''''^" + ^ ^ Gm(Q^) (38) 

where Q'^ = if — luq = —2M]s[UJq, and Ge{Q^) and Gm{Q^) stand for the electric and magnetic nucleon form factors 
respectively. The spin independent GPs p(^i^ii)o and p(*^i.A^i)o a,re directly proportional to the electric and magnetic 
polarizabilities of the nucleon respectively. A Fourier transform of their g-dependence allows us to map out the spatial 
distribution of the electric polarization and magnetization of the nucleon, as discussed in Ref. 

The first VCS experiment was performed at MAMI ,2,] and the response functions Plt and Pll - Ptt/s were 
extracted at = 0.33 GeV^. Going to higher Q^, the VCS process has also been measured at JLab and data have 
been obtained a.t = 1 GeV^ and = 1.9 GeV^ [^. To unambiguously extract the spin- independent GPs from 
experiment, requires to separate Pll from Ptt- This can be done by performing two unpolarized VCS experiments 
at a fixed value of q (or equivalently Q^) and by varying e (or equivalently the beam energy). Such experiments are 
planned at MAMI in the near future, making use of the beam energy upgrade. Besides extracting the electric CP 
p{Li,Li)o^ such an experiment will also allow for a measurement of the combination of the two spin GPs p(*fi.*fi)i 
and p(ii.A")i 

Until now, we discussed only unpolarized VCS observables. An unpolarized VCS experiment gives access to only 
three combinations of the six GPs, as given by Eqs. (|^ - (|2Hl- It was shown in Ref. 0| that VCS double polarization 
observables with polarized lepton beam and polarized target (or recoil) nucleon, will allow us to measure three more 
combinations of spin-flip GPs. Therefore a measurement of unpolarized VCS observables (at different values of e) 
and of three double-polarization observables will enable us to disentangle all six GPs. The VCS double polarization 
observables, which are denoted by AA4{h, i) for an electron of helicity are defined as the difference of the squared 
amplitudes for recoil (or target) proton spin orientation in the direction and opposite to the axis i {i = x, y, z), where 
the z-direction is chosen along the virtual photon momentum (see Ref. |0| for details). In a LEX, this polarized 
squared amplitude yields : 

^_^exp ^ ___2_ ^ ^ ^_^cxp ^ Q(^,) (3g) 

Analogously to the unpolarized squared amplitude of Eq. I|34l) . the threshold coefficients AAl^j', AA^^f are known 
due to the LET. It was found in Ref. 0| that the polarized squared amplitude AA^q'^'' can be expressed in terms 
of three new structure functions P£rp(g), Pjlp{q), and Pjjf,{q). These new structure functions are related to the spin 
GPs as |l[l3 : 

PLt = ^Gm(Q^)p("'^^)^ - ^^"g^,(Q2)p(*"'^")\ (40) 
P,^^ = _^QGm(Q^)p("'^^'^ + ^^^G^,(Q2)p(^"-^")\ (41) 
p'l^t = Y^GMiQ') (^P("'")i - y|^oP(^^^'"^^j . (42) 

While P£j. and P^^ can be accessed by in-plane kinematics (0 = 0°), the measurement of P^^ requires an out-of-plane 
experiment. 
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q=600MeV q'=111.5 MeV e=0.62 <t)=0" 





-150 -100 -50 



150 -100 -50 
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FIG. 2: VCS double-polarization asymmetry (polarized electron, recoil proton polarization along either the z- or x- directions 
in the cm. frame) in MAMI kinematics as function of the photon scattering angle. The dashed-dotted curves correspond to 
the BH+B contribution. The dotted curves are the lowest order HBChPT predictions from Refs. The dashed curves are 

obtained by including the next order HBChPT correction for all spin-flip GPs, as calculated in this work, while using the lowest 
order predictions for th e sp in-independent GPs. For comparison, we also show the results of a phenomenological dispersion 
relation calculation llB.ll8|| (solid curves). 



In Fig. 12 we show the results for the double polarization observables, with polarized electron and by measuring the 
recoil proton polarization either along the virtual photon direction (z-direction) or parallel to the reaction plane and 
perpendicular to the virtual photon (x-direction). These asymmetries directly depend on the response functions of 
Eqs. H4l)|) - (|42|l . One sees from Fig.[21that the double polarization asymmetries are quite large (due to a non-vanishing 
asymmetry for the BH -I- B mechanism). We furthermore compare in Fig. |2 the HBChPT predictions for the lowest 
order (non- vanishing) results for the spin GPs from Refs. 0,1^, with the next-order results for the spin GPs, as 
calculated in this work. For the spin independent GPs, which also enter in these asymmetries, we use in both cases 
the LO results, awaiting a NLO calculation for p(-E'i'i'i)o ^j^^j p(Mi,Afi)o^ 0-a& observes from Fig. [21 sizeable differences 
in the asymmetries between the leading order and next order HBChPT predictions for the spin GPs. For comparison, 
we also show in Fig.|21the results of the DR calculation of Refs. [l^[l8l |. Due to the smaller spin GPs in the dispersion 
calculation, the relative effect on the asymmetries is smaller in the DR approach. 

To quantify the relative effect on the asymmetries due to the spin GPs, we display in Fig.Othe double polarization 
asymmetries with the well known Bethe-Heitler -|- Born contribution subtracted. Another contribution which is well 
known is the anomaly contribution (due to a i-channel tt'^ pole) to the spin GPs, which is shown in Fig. |31 (lower panels). 
One furthermore sees from Fig. |21 (upper panels) that the effect of the spin GPs, excluding the anomaly contribution, 
on the double polarization asymmetries is of the order of 5 %. To distinguish between the different predictions, a 
double polarization VCS experiment is clearly called for. Although these double polarization observables are tough 
to measure, a first test experiment is already underway at MAMI |jT9j |. 
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q=600MeV q'=111.5 MeV £=0.62 ^=0" 
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FIG. 3: Deviation of the double-polarization VCS asymmetry from the BH+Born result, calculated within the LEX formalism in 
MAMI kinematics as function of the photon scattering angle. Upper panels: lowest order HBChPT predictions from Refs. 0,0 
(dotted curves); results including the next order HBChPT corrections for the spin-flip GPs, as calculated in this work, and 
using the leading order predictions for the spin-independent GPs (dashed curves); dispersion relation results 17^ (solid 
curves). The effect of the anomaly contribution (i.e. t-channel 7r° pole) is neglected in all the three model-predictions. In the 
lower panels, a comparison is shown between the dispersive predictions without the anomaly contribution (dashed curves) and 
with the anomaly contribution (solid curves). 



V. CONCLUSIONS 



In this work, we calculated the spin-dependent VCS amplitude in HBChPT at 0{p'^). In a low energy expansion of 
the outgoing photon energy, the spin-dependent VCS amplitude can be parametrized in terms of 7 generalized spin 
polarizabilities which characterize the response of the nucleon. Nucleon crossing symmetry combined with charge 
conjugation invariance leads to relations among these GPs, leaving four independent spin-flip GPs. We calculated 
all seven spin-flip GPs in HBChPT at 0{p'^) and checked the three crossing symmetry relations to this order. At 
order ©(p^), no unknown low-energy constants enter the theory, allowing us to make absolute predictions for all 
seven spin-flip GPs to this order. Three of the spin-flip GPs at next order have been extracted before in Ref. 6] by 
identifying the term linear in the outgoing photon energy in the LEX of the spin-dependent VCS amplitudes. In the 
present work, we calculated within HBChPT at ©(p*) the quadratic terms in the outgoing photon energy, in which 
the dipole GPs also appear. From these quadratic terms, we were able to extract the remaining four spin-flip GPs at 
next order in HBChPT. A noteworthy feature of our results is that we provide analytical expressions for all spin-flip 
GPs at next order, analogously as has been done at lowest order in Ref. |^. Another interesting feature arises from the 
crossing symmetry relations which connect the HBChPT amplitudes at different orders. As a result, our prediction 
for the GP p(^i'ii)i is formally a N^LO result, and could only be obtained from a HBChPT calculation at 0{p^), if 
one were not to use the crossing relations. 

Our predictions show sizeable corrections to the leading order (non-zero) results for the spin- flip GPs. The sign of 
these corrections is such as to move the leading order results in the direction of a phenomenological estimate of GPs 
based on dispersion relations. As the spin-flip GPs allow for absolute predictions at two successive orders in HBChPT 
without new low-energy constants entering, they may also be useful as a testing ground of chiral theories. It would 
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therefore be desirable to have direct experimental information on the spin- flip GPs to test the different predictions. 

One combination of the four spin-flip GPs can be extracted by doing an unpolarized VCS experiment at a fixed 
value of and by varying the beam energy. One can extract the other three spin-Hip GPs through VCS double 
polarization observables with polarized lepton beam and polarized recoil nucleon. We gave predictions of these double 
polarization asymmetries for experimental conditions accessible at MAMI. We found that the relative effect on these 
asymmetries due to the spin- flip GPs is of the order of 5 %. Although the measurement of these asymmetries requires 
dedicated experiments, such a measurement would challenge the parameter-free chiral predictions at next order for 
the spin-flip GPs presented in this work. 
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APPENDIX A: CROSSING SYMMETRY PROPERTIES OF THE VCS AMPLITUDES IN THE 

COVARIANT BASIS 

In this appendix, we study the properties of the amplitudes fi under nucleon crossing combined with charge 
conjugation, which transforms fi{q^,q ■ q',q ■ P) into fi(q^,q ■ q',—q ■ P)- First we expand the arguments of the 
amplitudes fi as : 



Mq',q-q\q-P) - J] C,, • (g^f • (g • q')" ' ' 



Now we obtain : 



g2 = 2AfjvcJo + 2t^o-^' + (-|^ + l)-cj'2 + C'(cj'3), 

Mat 

q-q' = 0+ {uio-q COS e) ■ uj' + 1-Lj'^ +0{uj'^), 

q-P = 0+{— +qcos0) -uj' + l-uj''^ + 0{uj'^). (Al) 



f^ = ^C:oo(2Afwa;o)", 

oc c« oo 

9^ = IZQio(2^^wWo)"^o + ^Qoi(2Mjv^o)"(2MAr - uj^) + ^Qoo ' a(2Mjv^o)""' • 2^o, 

a— a— a— 1 

oc oo 

= ^-C:io(2MAfC^o)"g + 5IOi(2^^w^o)"g, (A2) 

0=0 a=0 

where Qj,^ are the coefficients which start at rt = 0, and Qf,^, = 'Y^=q aT^ ' ^^abl ■ 

It was shown in Ref. Ji]] that nucleon crossing symmetry combined with charge conjugation leads to the following 
symmetry properties of the amplitudes : 

Mq^,q-q',q-P) = n{q^ , q ■ q' , -q ■ P), i = 1, 2, 5, 6, 7, 9, 11, 12, 

Mq^,q-q',q-P) = -Mq^ ,q ■ q\~q ■ P), z = 3,4,8,10. (A3) 
Furthermore, photon crossing leads to the symmetry relations at the real photon point : 

f,iO,q-q\q-P)^0, i = 7, 9. (A4) 
From these symmetry relations, we obtain for the coefficients entering the expansion of Eq. IjAljl : 

Coi = 0, 1 = 1,2,5,6,11,12, 

Cio = Coo = 0, z = 3,4,8,10, 

Qoi — ^oiQ — i = 7,9. 

(A5) 
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As a result, we obtain for the LEX of Eq. ifTH) : 

oo oo 

a— a— 1 

oo 

/l^ = ^-Qio(2M^'^o)"g, for 1 = 1,2,5,6,11,12, (A6) 



a=0 



furthermore 



h = 0, = ^Coi(2Mjvc^o)'^(2AfA,-c^o), 

oo 

= IICoi(2A/wwo)"g, for z = 3,4,8,10, (A7) 



a=0 



and 



5» = ^Qio(2Mivwo)''wo + ^Qoo • a(2Mjvwo)""' • 2wo, 

a— 1 a— 1 

oo 

/iz = ^-Qio(2MjvWo)"9, for i = 7,9. (A8) 

a=l 

The above results in turn yield relations for the coefficients in the heavy baryon expansion as in Eqs. H2()I21|I . Since 
in the heavy baryon expansion, ujq is counted as n = 1 and 2Mj^uJo as n = 0, we obtain : 

.9^ = 9i"' - .9^ = .9^"^ - - af' = = = 0. (A9) 

Moreover we know 

9f(0) , j-(O) _ .(0) _ „ 

^ie + ^9 ~ Ji2 — 

because otherwise the amplitudes Ag and A12 would start at n = — 1 which is impossible. Therefore 

2/i°^ + h'' - /!? = 0, (AlO) 

and analogously 

2/^(°) + 4°) = h^^ = 0. (All) 

From = it follows that C^io"""' = 0, and from j'^ = it follows that C^oo"~°'' — 0- Therefore we obtain 
.912'' ~ 0- Similarly, one can show that : 

/if = /if = /jf") = hf^ = 0. (A12) 
otherwise gs, 94, gs^ 9io will start at n = — 1. Lastly, Eq. HA7|I yields the relations : 

/ir^=ffr-^, for * = 3,4,8,10. (A13) 
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